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7.  Background/Scope  of  Effort:  The  physics  of  interaction  between  key  materials  and 
electromagnetic  radiation  at  terahertz  (THz)  frequencies  of  1012  to  1013  Hz  has  been  considered 
promising,  on  two  counts,  for  the  detection  of  concealed  explosives.  First,  characteristic 
absorption  features  in  reflection  spectra  of  explosives  may  provide  a  basis  for  discrimination  from 
other  materials.  Second,  common  obscuring  materials  (e.g.,  cellulose,  plastic)  are  translucent  to 
THz  radiation,  so  detection  and  imaging  of  spectral  signatures  inside  containers  have 
considerable  potential.  However,  realization  of  this  potential  is  both  mediated  by,  and 
complicated  by,  effects  of  scattering.  Scattering  from  the  granular  structures  and  rough  surfaces 
of  actual  explosives  may  distort  or  obscure  the  spectral  features.  This  is  a  general  issue  with  all 
THz  techniques  for  imaging  and  detection  of  a  strongly  scattering  target.  It  is  crucial  to 
understand  the  volume  and  rough  surface  scattering  processes,  to  understand  the  interaction 
between  scattering  and  spectral  absorption,  and  to  develop  new  theoretical  models  and  new 
technology  to  extract  the  spectral  absorption  information  in  the  presence  of  strong  scattering.  On 
the  other  hand,  scattering  may  make  spectral  signature  more  easily  detectable,  because  diffuse 
scattering  can  send  the  scattered  energy  over  a  wide  range  of  backscattering  angles,  much  more 
convenient  for  detection  then  the  case  of  specular  reflection  from  a  planar  target,  where  the 
geometry  would  need  to  be  tightly  controlled  to  permit  detection  of  spectral  signatures. 

8.  Summary/ Abstract:  The  research  conducted  under  this  project  was  a  joint  effort 

between  Applied  Physics  Laboratory,  University  of  Washington  (APL-UW),  and  Portland  State 
University  (PSU).  The  focus  at  APL-UW  was  on  1)  developing  an  experimental  setup  for 
variable  source  and  detector  angles  and  variable  linear  polarization,  and  then  implementing  it  on  a 
commercial  THz  spectroscopic  system  at  PSU,  2)  expanding  the  current  THz  time  domain 
spectroscopic  system  at  APL-UW  to  cover  broader  bandwidth  by  developing  an  air  plasma  THz 
source,  and  3)  developing  analysis  techniques  that  allow  spectroscopic  signatures  to  be  recovered 
in  the  presence  of  appreciable  background  noise,  a  situation  that  can  be  expected  when  volume 
and  rough  surface  scattering  processes  are  important.  The  long-term  goal  of  this  project  has  been 
to  lay  the  groundwork  for  an  integrated  spectral  and  imaging  methodology  on  which  to  base 
technology  development. 

9.  Technical  Contents  and  Accomplishments 

9.1.  Implementation  of  an  automated  measurement  system 

The  task  of  the  first  program  year  was  to  implement  automated  variable  transmitter  and 
receiver  angles  and  variable  linear  polarization  measurements  on  a  commercial  THz 
spectroscopic  system  at  Prof.  Lisa  Zurk’s  laboratory  at  Portland  State  University.  Prof.  Zurk  was 
a  co-PI  of  our  previous  project  in  the  SAAET  program,  and  although  she  now  leads  a  separate 
SAAET  project,  her  project  and  our  project  have  been  collaborative  in  nature  with  similar  goals. 

When  a  THz  wave  is  scattered  on  a  rough  surface,  the  incident  THz  wave  is  re-distributed  over 
a  wide  range  of  angles  and  the  state  of  polarization  of  the  scattered  wave  can  be  different  from 
that  of  the  incident  wave.  It  is  necessary  to  perform  THz  scattering  measurements  at  various 
angles  and  polarizations  to  accurately  describe  the  scattering  properties  of  a  test  sample.  PSU 
acquired  a  commercial  THz  spectroscopy  system  made  by  Picometrix.  It  is  a  fiber  laser  based 
system  that  offers  great  flexibility  and  allows  the  transmitter  and  receiver  heads  to  be  moved 
freely  and  be  arranged  with  any  angles  of  incidence  and  reflection  with  respect  to  the  test  sample. 
Because  scattering  is  a  stochastic  process  in  a  random  medium  or  at  a  randomly  rough  surface,  a 
large  amount  of  data  is  needed  to  develop  and  validate  statistical  models.  To  accomplish  this 
requires  a  system  of  computer  controlled  motion  stages  in  precisely  coordinated  movement. 
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The  system  of  motion  stages  was  designed  and  constructed  at  APL-UW  (Figure  1).  It  consists 
of  a  stack  of  two  rotation  stages  to  rotate  the  sample  and  the  receiver  (the  transmitter  is  stationary 
and  the  change  of  angle  of  incidence  is  achieved  by  rotating  the  sample),  an  X-Y  translation  stage 
for  in-plane  motion  of  the  sample  for  different  realizations  of  the  scattering  measurement,  and 
two  rotation  stages  to  rotate  the  transmitter  and  receiver  heads,  respectively.  After  it  was 
assembled  and  tested  it  was  transferred  to  PSU  and  installed  with  the  Picometrix  transmitter  and 
receiver  heads  to  form  an  automated  THz  scattering  measurement  system.  The  APL-UW  team 
trained  the  PSU  team  on  how  to  interface  the  system  with  computers  and  LabView  software,  and 
how  to  operate  the  system.  The  APL  team  also  provided  technical  support  throughout  the 
duration  of  the  project. 


Figure  1.  A  view  of  the  automated  motion  stages  installed  on  the  PSU  THz  spectroscopy  system. 


The  automation  of  the  measurement  system  enabled  the  two  groups  to  do  extensive  and 
systematic  studies  of  THz  scattering.  Without  the  automation  system,  setting  up  different 
measurement  configurations  and  collecting  large  amounts  of  data  had  to  be  done  manually  and 
was  prone  to  human  error.  The  system  proved  to  be  highly  useful  for  research  during  the 
subsequent  program  years. 

9.2.  Development  broadband  air  plasma  THz  source 

The  approach  taken  to  develop  a  broader  band  THz  source  was  to  utilize  an  air  plasma  for  the 

source.  The  nonlinear  process  by  which  the  third  order  optical  susceptibility,  ^(3) ,  contributes  to 
the  generation  of  THz  radiation  is  known  as  four-wave  mixing.  It  was  first  shown  by  Cook  and 
Hochstrasser  [1]  in  2000  that  the  air  plasma  generated  by  ionization  of  air  molecules,  when  a 
femsto-second  laser  beam  is  focused  in  the  ambient  air,  exhibits  a  very  strong  third  order 
susceptibility  tensor,  which  can  be  exploited  to  efficiently  generate  intense  coherent  THz 
radiation.  The  process  uses  the  nonlinear  interactions  of  the  fundamental  wavelength  photon, 
which  also  produces  the  air  plasma,  with  its  second  harmonic  photons  generated  in  a  BBO  crystal. 
We  have  implemented  a  four- wave  mixing  THz  source  in  our  THz-TDS  setup.  Figure  2  shows 
photographs  of  the  air  plasma  size  and  intensity  generated  using  our  amplified  femtosecond  laser 
beam  with  average  power  of  800  mW,  center  wavelength  of  790  nm,  and  pulse  length  of  50  nm  or 
45  fs,  produced  by  focusing  through  lenses  with  various  focal  lengths.  Note  that  the  camera 
exposure  time  has  been  significantly  increased  from  the  strongest  lens  to  the  case  with  the  longest 
focal  length  in  order  to  resolve  and  fully  capture  the  ionized  region  in  the  ambient  air. 
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Fig.ure  2.  Effect  of  varying  the  focal  length  f  on  the  size  and  intensity  of  the  air  plasma  produced, 
using  the  femtosecond  amplified  laser  system  at  UW-APL. 

It  has  been  shown  experimentally  that  tighter  focusing  produces  a  more  intense  air  plasma; 
however,  for  physical  considerations  such  as  the  phase  matching  conditions,  the  length  of 
interaction  between  the  fundamental  and  its  second  harmonic,  the  profile  of  the  air  plasma  region, 
as  well  as  the  convergence  angle  of  the  THz  radiation,  longer  focusing  results  in  stronger 
terahertz  emission  [2-6].  Figure  3  shows  the  higher  power  and  broader  bandwidth  of  the  terahertz 
waves  generated  using  this  technique  in  a  N2  purged  air  plasma  source  vs.  optical  rectification  in 
ZnTe  for  comparison.  We  have  also  explored  using  various  EO  crystals  to  extend  the  detection  of 
bandwidth  of  our  system  to  support  future  experiments.  Figure  3  also  shows  that  a  broader 
bandwidth  can  be  achieved  using  a  GaP  detector  compared  to  ZnTe  crystals.  This  experimental 
capability  would  allow  measurements  of  incoherent  THz  spectroscopy  from  rough  surface 
substances  to  be  expanded  to  much  higher  THz  frequencies. 

9.3.  Implementation  of  multi  resolution  analysis  for  identifying  spectroscopic  features 

When  attempting  to  identify  spectroscopic  features  in  a  THz  signal  scattered  from  a  target 
with  a  rough  surface  or  with  volume  heterogeneity,  the  signal  is  likely  to  be  embedded  in 
background  noise  produced  in  the  scattering  process  [7].  This  can  be  illustrated  by  observing  how 
the  strength  of  a  spectral  feature  decreases  as  the  roughness  on  a  test  sample  increases.  We  chose 
as  our  test  material  a-lactose  monohydrate,  whose  absorption  signatures  at  0.54,  1.20,  and  1.38 
THz  are  well  known.  Controlled  levels  of  surface  roughness  were  inscribed  in  the  pellets  from 
three  different  grades  of  sandpaper  sheets,  designated  as  400,  150,  and  PI 80  grit.  We  also 
measured  the  roughness  parameters  of  the  sandpaper  independently  to  support  theoretical 
modeling.  We  then  used  our  THz  time-domain  spectrometer  (optical  rectification  in  ZnTe)  to 
measure  scattered  radiation  in  the  specular  direction  for  separate  patches  on  the  rough  surface, 
from  which  we  computed  the  mean  power,  individual  incoherent  scattered  powers  (power  for  a 
single  realization  minus  the  mean  coherent  power),  and  the  mean  incoherent  scattered  power,  all 
as  functions  of  THz  frequency.  A  sample  of  results  is  shown  in  Figure  4,  where  the  roughness 
increases  as  the  grit  size  decreases.  Figure  4  shows  the  Fourier  amplitude  of  the  reflected  THz 
field  as  a  function  of  frequency,  and  the  lactose  spectroscopic  feature  of  primary  interest  is  at  0.54 
THz. 

Important  spectroscopic  features  can  be  made  more  easily  visible  by  plotting  the  first 
derivative  of  the  spectral  amplitude,  as  shown  in  Figure  5  where  the  frequency  range  as  been 


Final  Report  for  ONR  CODE  30  SAAET  6.1  Basic  Research  Effort 


reduced  and  centered  on  the  spectral  feature  of  interest.  Note  that  as  the  roughness  increases,  the 
visibility  of  the  absorption  feature  decreases  and  is  effectively  lost  for  the  two  roughest  cases.  We 
also  found  the  frequency  dependence  of  coherent  power  to  be  in  close  agreement  with  simple 
theoretical  modeling  based  on  the  Kirchhoff  approximation  (inset). This  information  was  used  to 
detrend  the  spectral  amplitudes  prior  to  the  wavelet  analysis  to  be  described  shortly. 


Figure  3.  A  comparison  between  the  spectral  content  of  the  terahertz  field  generated  using  both 
the  four-wave  mixing  process  in  an  air  plasma  and  optical  rectification  in  a  ZnTe  crystal.  The 
terahertz  fields  are  detected  by  the  EO-sampling  method  using  both  a  ZnTe  and  GaP  crystal. 


10-3i - 1 - 1 - 1 - 1 - 1 
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Figure  4.  Fourier  amplitude  of  the  reflected  THz  field  for  smooth  and  rough  lactose  pellets. 
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The  result  in  Figure  5  illustrates  the  need  for  signal  processing  techniques  that  can  detect  the 
presence  of  spectral  features  that  might  be  obscured  by  fluctuations  in  the  background  level 
(which  we  refer  to  as  noise).  Motivated  by  this  need,  we  have  developed  a  method  using  wavelet 
decomposition  of  the  frequency-domain  signal  that  successfully  recovers  the  location  of  the 
spectral  feature  (resonance)  even  in  the  roughest  case  shown  in  Figure  5  (P80  grit).  Wavelet 
methods  are  particularly  suitable  for  analysis  of  data  series  with  localized  features,  whose 
locations  are  not  known  a  priori ,  as  opposed  to  Fourier-based  spectral  analysis  techniques  that 
work  best  with  stationary  functions.  The  wavelet  method  works  essentially  by  separating 
variations  in  the  frequency-domain  signal  into  a  sequence  of  scale  sizes  (in  this  case  frequency 
scales)  with  the  effect  of  separating  out  the  fluctuations  in  the  background  caused  by  incoherent 
scattering  from  the  specific  scale  of  the  resonance,  thus  providing  a  statistically  sound  basis  for 
detection  of  the  resonance  even  when  it  is  obscured  by  incoherent  scattering  effects. 


v>  x  10"'2 


Figure  5.  First  derivative  of  reflection  spectral  amplitude  of  lactose  pellets  from  Figure  4. 

A  detailed  discussion  of  wavelet  methods  in  general  and  the  particular  choices  made  for 
application  to  analyzing  terahertz  reflection  spectra  would  require  an  extensive  treatment,  even  to 
clarify  the  technical  terms  used  in  describing  the  analysis.  Such  an  extensive  treatment  will  not  be 
included  here  since  the  details  are  available  in  the  literature.  A  very  useful  treatment  on  wavelet 
methods  in  general  can  be  found  in  the  monograph  by  Percival  and  Walden  [8].  Details  on  our 
analysis  are  also  available  [9-11]. 

An  illustration  of  how  the  wavelet  method  can  separate  the  variations  in  the  frequency-domain 
signal  into  a  sequence  of  frequency  scale  sizes  is  shown  in  Figure  6,  where  the  original  detrended 
spectral  data  series  from  the  roughest  case  in  Figures  4  and  5  (roughness  obtained  with  the  P80 
grit)  is  shown  at  the  top  and  labeled  X.  The  vertically  offset  plots  below  result  from  a  wavelet 
decomposition  of  data  series  X  into  a  sequence  of  components  ranging  from  the  largest  frequency 
scale  size  at  the  top  to  the  smallest  at  the  bottom.  When  the  eight  offset  plots  Dx  to  D8  are  added 
together  they  reconstruct  the  original  data  series  X.  The  absorption  feature  of  lactose  near  0.54 
THz  is  evident  in  D2  and  D3 . 

In  Figure  6  all  of  the  offset  plots  have  the  same  vertical  scale,  causing  the  signature  of  the 
absorption  feature  at  0.54  THz  to  become  small  for  the  shortest  scale  sizes,  though  the  signal  to 
noise  ratio  remains  good.  Thus,  expanding  the  vertical  scale  for  the  shortest  scale  sizes  is 
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beneficial.  This  is  shown  in  Figure  7  for  data  obtained  with  the  air  plasma  source  illuminating 
roughened  lactose  (P80  grit). 
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Figure  6.  Separation  of  the  original  rough  lactose  (P80  grit)  detrended  spectral  data  series  X 
into  a  sequence  of  spectral  scale  sizes  using  the  wavelet  method.  Vertically  offset  plots  show 
components  of  the  original  data  series  ranging  from  large  to  small  scale  size  (top  to  bottom).  The 
absorption  resonance  of  lactose  near  0.54  THz  is  evident  in  D2  and  D3 .  The  vertical  axis  is 
dimensionless. 
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Figure  7.  Separation  of  the  original  rough  lactose  (P80  grit)  detrended  spectral  data  series  X 
into  a  sequence  of  spectral  scale  sizes  using  the  wavelet  method.  The  air  plasma  source  was  used 
for  this  case,  and  only  the  four  smallest  scale  components  are  shown. 
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The  data  series  components  shown  in  Figure  7  correspond  to  the  four  lowest  shown  in  Figure 
6,  though  the  notation  labeling  the  components  has  been  changed.  It  is  evident  that  for  shortest 
spectral  amplitude  scale  (Dx  in  Figure  6,  TWi  in  Figure  7)  the  absorption  resonance  feature  of 
lactose  near  0.54  THz  has  the  highest  signal  to  noise,  which  would  facilitate  an  automated 
detection  algorithm.  Additional  absorption  features  also  appear  in  the  region  above  1  THz  for  this 
air  plasma  source  generated  data  series. 

Another  issue  that  would  need  to  be  addressed  for  an  automated  detection  algorithm  is  the 
detrending  step  that  was  used  to  remove  the  main  decrease  in  signal  level  as  the  frequency 
increases,  shown  in  Figure  4.  The  detrending  operation  in  that  case  was  performed  based  on  the 
frequency  dependence  of  the  coherent  reflection  loss  from  a  rough  surface,  and  then  numerically 
fitting  that  trend  to  the  data  as  shown  in  the  inset  in  Figure  5.  It  would  be  much  better,  however, 
to  have  a  universal  detrending  algorithm  that  does  not  depend  on  fitting  a  theoretical  scattering 
prediction  to  the  data.  It  happens  that  the  wavelet  approach  being  used  imposes  circular  boundary 
conditions  on  the  data  series,  which  is  one  reason  that  an  initial  detrending  step  is  made  before  the 
wavelet-based  decomposition  is  performed.  An  approach  was  developed  to  avoid  the  detrending 
step  altogether  by  using  the  original  data  series  as  the  first  quarter  of  an  artificial  periodic 
extension  such  that  the  resultant  extended  data  series  would  repeat  itself  by  the  completion  of  the 
full  period.  Figure  8  illustrates  this  procedure. 
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Figure  8.  (a)  A  periodic  extension  of  the  P80  lactose  spectral  amplitude  (solid  line)  along  with 
the  so-called  smooth  part  (dotted  line)  which  carries  the  dominant  large-scale  structure,  (b)  The 
wavelet  decomposition  of  the  data  series,  similar  to  Figure  6,  after  removal  of  the  smooth  part. 
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Figure  8(a)  shows  the  periodic  extension  of  the  original  P80  grit  spectral  amplitude  along  with  a 
“smooth  part”  that  also  is  obtained  as  part  of  the  wavelet  processing  scheme.  After  its  removal, 
the  decomposition  into  a  series  of  spectral  scale  sizes  proceeds  _as  before,  and  the  spectral 
absorption  feature  along  with  its  periodic  extension  is  visible  in  the  D2  and  D3  components.  If  the 
vertical  axis  were  expanded  as  in  Figure  7,  then  the  absorption  feature  would  be  just  as  visible 
in  Di  as  well.  Following  this  procedure  an  automated  algorithm  could  easily  by  implemented  for 
detecting  specific  spectral  features. 

Another  issue  that  is  worth  understanding  is  the  extent  to  which  the  wavelet  analysis  can 
distinguish  closely  spaced  spectral  features.  To  study  this  topic,  a  realistic  model  for  an 
absorption  feature  was  added  to  a  measured  spectral  amplitude  obtained  from  smooth  lactose.  The 
result  is  a  spectral  amplitude  with  two  absorption  features,  one  at  the  native  0.54  THz  of  a- 
lactose,  the  second  at  a  frequency  that  was  varied  in  steps  from  0.5  to  0.6  THz.  The  results  are 
shown  in  Figures  9-13.  As  with  Figure  7,  only  the  four  smallest-scale  components  obtained  from 
the  wavelet  analysis  are  shown,  and  the  components  are  also  labeled  as  in  that  figure.  It  can  be 
seen  that  as  the  second  peak  moves  from  0.5  to  0.6  THz,  except  for  when  it  is  located  at  0.55  THz, 
the  existence  of  two  separate  features  can  be  clearly  identified  When  the  twin  feature  is  placed  at 
0.55  THz  in  Figure  11,  very  close  to  the  native  one  at  0.54  THz,  it  is  still  possible  to  observe  that 
the  structure  is  more  complex  than  obtained  from  a  single  absorption  resonance.  Therefore,  the 
possibility  remains  that  with  suitable  modeling  input,  an  automatic  algorithm  could  be  produced 
that  would  recognize  this  special  situation  and  flag  it  for  more  detailed  analysis. 


Frequency  (THz) 

Figure  9.  The  THz  reflection  spectra  of  an  artificial  material  with  the  additional  simulated  twin 
absorption  feature  of  lactose  at  0.50  THz. 


In  summary,  the  wavelet  method  that  has  been  developed  shows  considerable  promise  in 
identifying  spectral  features  that  may  be  embedded  in  background  noise.  A  patent  application  has 
been  filed  through  the  University  of  Washington  Office  of  Technology  Transfer  on  this  method, 
which  likely  has  application  well  beyond  THz  reflection  spectroscopy. 
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Figure  10.  The  THz  reflection  spectra  of  an  artificial  material  with  the  additional  simulated  twin 
absorption  feature  of  lactose  at  0.52  THz. 


Figure  11.  The  THz  reflection  spectra  of  an  artificial  material  with  the  additional  simulated  twin 
absorption  feature  of  lactose  at  0.55  THz. 
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Frequency  (THz) 


Figure  12.  The  THz  reflection  spectra  of  an  artificial  material  with  the  additional  simulated  twin 
absorption  feature  of  lactose  at  0.58  THz. 
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Figure  13.  The  THz  reflection  spectra  of  an  artificial  material  with  the  additional  simulated  twin 
absorption  feature  of  lactose  at  0. 60  THz. 
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